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Measuring parameters of AGN central engines with very 
high energy 7-ray flares 
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ABSTRACT 

We discuss a " compact source" model of very high energy (VHE) emission from blazars 
in which the variability time is determined by the blazar central engine. In this model 
electron or proton acceleration close to the supermassive black hole is followed by 
the development of electromagnetic cascade in a radiatively inefficient accretion flow. 
Assuming such a model for the TeV blazar PKS 2155-304, we show that the variability 
properties of the TeV 7-ray signal observed during a bright flare from this source, 
such as the minimal variability time scale and the recurrence period of the sub- flares, 
constrain the mass and the angular momentum of the supermassive black hole. 



Key words: gamma-rays: theory, galaxies: nuclei, radiation mechanisms: non- 
thermal, black hole physics, BL Lacertae objects: individual: PKS 2155-304 



1 INTRODUCTION 

Recent observation of fast variabil ity of TeV 7-ray emis- 
sion from several TeV blazars (jAharonian et al . 2007; 



lAlbert et all l2007h challenges the conventional model in 
which the TeV 7-rays are supposed to be produced at large 
distances from the blazar central engine, the supermassive 
black hole. Within this conventional model the 7-ray emit- 
ting blobs are assumed to travel to the parsec-scale dis- 
tances along the AGN jet before emitting in the TeV en- 
ergy band; it is believed that otherwise the 7-ray emission 
would be strongly absorb ed in the accretion flow (see e.g. 
iGhisellini fc Madaul l|l99rj )). Since the radiative cooling time 
of the TeV emitting electrons is typically shorter than the 
time of propagation from the central engine to the TeV emis- 
sion region, it is usually assumed that the TeV emitting 
particles are produced via shock acceleration locally in the 
emission region, rather than in the AGN central engine. 

At the same time, the observed short variability time 
scales i var ~ a few minutes indicate that the TeV emission 
comes from very compact regions having the size in the co- 
moving frame Ax' < 8(1 + z) -1 c£ V ar, where <5 is the bulk 
Doppler factor and z is the source redshift. This implies that 
in the static frame the longitudinal size of the TeV emitting 
region is 



At' 

Ax=^<(1 + 2 )- 
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tv 
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(1) 



where we assume that the bulk Lorentz factor V ~ 8. This 
is comparable to the minimal possible scale set up by the 
gravitational radius of the central supermassive black hole 

R g = GMbh/c 2 ~ 1.5 x 10 12 [M BH /1O 7 M ] cm . (2) 

Even if one assumes that the TeV emitting plasma blobs are 
produced close to the black hole, which would explain their 
initially small size, it is not clear how the blobs propagat- 
ing downstream the relativistic jet can retain this size up 
to large distances, unless they have unreasonably large bulk 
Lore ntz factors. This probl em has recently led to a sugges- 
tion (Bcac lman et al.ll2008f ) that the TeV flares may be not 
triggered by the black hole but rather are results of enhanced 
emission intrinsic to the jet. 

However even in that case, to explain the observed rapid 
variability, the TeV 7-ray emitting blobs have to travel with 
quite large bulk Lorentz factors, r> 1. Such Lorentz factors 
are in contradiction with the radio observations of the mo- 
tion of hot spots in the parsec-scale jets. Moderate apparent 
speeds of the blobs of the parsec-scale jets, revealed by radio 
observations, combined with an estimate of the number of 
parent objects of TeV blazars, would give much smaller va l- 
ues of the bulk Lorentz factors, r ~ 1 ijHenri &: Sa uec 2006). 
For example, to explain the fast variability of the July 2006 
TeV flare of PKS 2155-304 jAharonian et alj|2007t ) the bulk 
Lorentz factor required by the mechanism of TeV emission 
in the parsec-scale jet shou ld be T > 50 l|Aharonian et all 
2007; Bcg efman et al.ll200Sl) . At the same time, the direct 
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observations of the apparent velocity of hot spots in PKS 
2155-304 jet at the projected distance (l-r2) parsecs give 
the value u app = (0.9 ± 0.3)c (jPiner et al.ll2008h . Assuming 
that the viewing angle is not too small, 9 > 1°, this yields 
r < 10 at the distance of a few tens of parsecs from the 
central engine. 

Both the problems of the fast variability and of the 
small observed Lorentz factors at parsec distances could be 
naturally resolved if the site of the VHE 7-ray production 
is located closer to the AGN central engine, at significantly 
sub-parsec distances. If the VHE emitting region is moving 
relativistically toward the observer with a bulk Lorenz factor 
r, the variability time scale limits the distance R of the 7-ray 
produ ction site from the central engine (see e.g. lCelotti et al] 

Gil)), 
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An immediate difficulty is, however, that at such distances 
the accretion fl ow onto the black hole can be opaque to the 
7-rays (see e.g. iBlandford fc Levinsonl ljl995h ). 

The problem of opacity of the compact source does not 
arise in the case of low-luminos ity AGNs that accr ete at sig- 
nificantly sub-Eddington rates ijCelotti et al.lll998l ). In these 
sources the accretion flow is described within the framework 
of t he radiatively in e fficient accretion flow (RIAF) mod- 
els ([Rees et all 1 19821 : iNaravan fc Yil 1 1994 Il995l : iNaravanl 
2002) in which most of the gravitational energy extracted 
from the accreted matter is converted into internal en- 
ergy, rather than into radiation. The possibility of escape 
of the VHE 7-rays from the vicinity of the AGN central en- 
gine is best illustrated by the nearby low-luminosity radio 
galaxy M87, which was recent ly found to be a source of the 
variable TeV 7-ra y emission {Ah aronian et al ] |2003l . 120061 : 
lAlbert et al. 20081). most probably coming from a compact 
source I Neronov fc Aharonianl 120071 : IRieger fc Aharonianl 
l2008bl ). 

In the compact source model the VHE 7-ray emission 
is triggered by high-energy particles accelerated close to 
the black hol e via o ne of the possib l e mec h anisms (see 
e.g. lLovelacd lll97rj): ILesch fc Pohll dl992T): iKardashev 



1 1995Ti : iBednarek fc Protherod J 1999)1: INeronov et al 
1 20021): iNeronov fc Semikozl (I2002T): iNeronov et all l|2005h : 



IRieger fc Aharonianl (|2008al ): INeronov et all i|2008D ). In 
this case the spectral and timing characteristics of the 
VHE emission are directly linked to the physics of the 
processes taking place close to the supermassive black hole, 
which naturally explains the variability of the signal on the 
shortest possible time scale. 

Within the AGN unification scheme, the TeV blazars 
(high-energy peaked BL Lacs) are assumed to be the beamed 
versions of the low- luminosity radio galaxi es similar to M87 
jBrownd fl983l : iGiroletti et alj liooj . I2OO6I ). Since the only 
difference between the TeV blazars and the low-luminosity 
radio galaxies is their orientation with respect to the line of 
sight, the compact source model can be applicable also in 
the case of TeV blazars. 

In what follows we adopt this point of view and de- 
velop a compact source model of high-energy activity of TeV 
blazars. We demonstrate that within such compact source 
model the characteristics of the fast-variable VHE emission 
can be used to constrain the parameters of the AGN cen- 



tral engine, in particular, the density and luminosity of the 
accretion flow, the black hole mass and spin. We illustrate 
this possibility on the particular example of the bright PKS 
2155-304 flare detected by the HESS telescope in July 2006 
|Aharonian et af]|2007i ). This flare consists of a number of 
well-pronounced sub-flares which exhibit quasi-periodic re- 
currence. We show that the rise time and the recurrence 
period of the sub-flares can be directly related to the light- 
crossing time and to the period of rotation over the last 
stable orbit around a Mbh ~ 10 7 Mq black hol^E 

The paper is organized as follows. In Sec. [5] we discuss 
the qualitative features of the model, including particle ac- 
celeration and propagation through the RIAF environment. 
The possibility of a new interpretation of the observational 
data within such a model is demonstrated in Sec.|3]where we 
find the constraints on the black hole mass and angular mo- 
mentum imposed by the timing analysis of the bright TeV 
flare of PKS 2155-304. We summarize our results in Sec. [4] 



2 COMPACT SOURCE MODEL OF THE TEV 
BLAZARS 

In the compact source model the VHE 7-ray emission region 
is assumed to be situated in the vicinity of the AGN cen- 
tral engine, rather than at parsec-scale distances. This fact 
implies two main differences of this model from a generic 
" relativistically moving blob" model of VHE emission from 
blazars. First, the high-energy particles responsible for the 
VHE emission can be injected into the 7-ray emission re- 
gion by the AGN central engine, rather than only by a 
process (shock acceleration) intrinsic to the blob. Second, 
the characteristics of the 7-ray emission are determined not 
only by the intrinsic properties of the blob, but also by the 
effects of propagation of the high-energy particles through 
the matter and radiation environment created by the ac- 
cretion flow. However, as we discuss below, the compact 
source model to some extent includes the blob model: the 
electromagnetic cascade, which develops as a result of the 
propagation of high-energy particles through the accretion 
flow environment, leads to the formation of a relativistically 
moving blob of secondary cascade particles. In the following 
sub-sections we consider general features of particle acceler- 
ation and propagation in the compact source model of TeV 
blazars. 



2.1 Particle acceleration and "direct" 7-ray 
emission from the acceleration region 

A number of mechanisms of particle acceleration in the 
vicinity of the central bl ack hole have b een proposed 
in the literature ( s ee e.g . lLovelacd (Il976h: ILesch fc Pohll 



( 19921): 
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Kardashev (11995 
d2QQ2|); 



Neronov et al 



Neronov et al 




.Bednarek fc Protherod |l99L. 
INeronov fc Semikozl (12003) 
IRieger fc Aharonianl (|2008al ) 
In this sub-section we summarize 
some common features of the acceleration models which 



1 This mass estimate is di fferent from the value ~ IO^Mq quoted 
bv lAharonian et al 1 ll2007ly We will comment on this discrepancy 
in Sec. 13.31 
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are related to the fact that the acceleration proceeds in 
a compact region of the size comparable to the black 
hole horizon and are independent of the details of the 
acceleration mechanism. 

A conventional dimensional estimate of maximal possi- 
ble energies of particles of charge e accelerated in a region 
of the size R ~ R g with the magnetic field B is 



B 



10 4 G 



Mb 



10 7 M ( 



eV 



(4) 



where k ^ 1 is the efficiency of a particular accelera- 
tion mechanism. This maximal energy is, however, not al- 
ways achieved because of the strong energy losses expe- 
rienced by the accelerated particles. The "minimal" en- 
ergy loss channel is the loss on the curvature radiation 
llLevinsonll2"000l ; lAharonian fc Neronovll2005l ; | Neronov et al.l 
2008), which limits the particle energies to 
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(5) 

where m is the particle mass. Here we normalize the mag- 
netic field to the value 10 4 G typical for the central engine of 
an AGN with a 10 7 Mq black hole. Additional energy losses, 
caused by the interactions of the accelerated particles with 
the matter and radiation backgrounds produced by the ac- 
cretion flow (see Sec. l2.2l below). lead to further reduction of 
the maximal attainable energy. If the magnetic field in the 
acceleration region is not ordered, the energies of the acceler- 
ated particles are reduced because of the strong synchrotron 
energy loss. 

Radiative energy losses which accompany particle ac- 
celeration (curvature, synchrotron, inverse Compton) result 
in the 7-ray emission directly from the acceleration region. 
If the resulting 7-rays are not completely absorbed during 
their propagation through the photon background created 
by the accretion flow, this "direct" 7-ray emission can, in 
principle, provide an observable signature of the compact 
source model. We will return to this point in Sec. 12.41 

Since the maximal energies of particles are determined 
by the balance of the acceleration and energy loss rates, all 
the work done by the electric field is dissipated via the avail- 
able energy loss channels. For each charged particle the en- 
ergy loss rate is about the acceleration rate, dE/dt ~ enBc, 
where we assume that the accelerating electric field strength 
is ~ kB. The maximal possible density of particles in the 
acceleration region, n q ~ nB/eR, is determined by the con- 
dition that a charge redistribution cannot screen the electric 
field. One estimates the total power extracted from the ac- 
celeration region by multiplying the energy loss rate of each 
particle on the maximal possible particle density and on the 
volume of the acceleration region, V ~ Rt\ this yields 



n q R g (dE/dt) 



(6) 
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Comparing this estimate to the apparent luminosity of a 
bright TeV flare (found under the assumption that the ra- 
diaiton is isotropic), Li so ~ 10 46 erg/s, one concludes that 
the energetics of such a flare is compatible with the com- 
pact source model if the TeV emission is beamed into the 
solid angle Q/4ir < 10 -3 . The emission from the compact 



source is expected to be anisotropic, with the direction of 
the beam of high-energy particles and/or photons set up by 
the direction of the magnetic field. 

A fraction of the work done by the accelerating electric 
field is carried away by the flux of accelerated particles (a 
relativistic wind). In the case of electrons, this fraction is 
usually negligible. In the case of protons or heavy nuclei the 
power of the relativistic particle wind, P p , can be estimated 
as 



Pr, ~ n q R 2 g c- min{_E cur , _B max } 

->42 



10 



(7) 



This order-of-magnitude estimate does not take into account 
a possible special geometry of the accelerating field. If, for 
instance, particles are accelerated in parallel magnetic and 
electric fields in the polar cap regions of black hole mag- 
netospheres, the curvature radii of particle trajectories can 
be somewhat larger than R g and, respectively, the maximal 
attainable energies and the pow er of the particle beam can 
be higher (jNeronov et alJl2008h . 

The high-energy particles and radiation emitted from 
the acceleration region are injected into the accretion flow. 
The form of the 7-ray signal produced by the interaction 
of the accelerated particles with the radiation coming from 
the accretion flow depends sensitively on the geometry of 
the central engine. The precise determination of this signal 
should involve detailed numerical modeling which is beyond 
the scope of the present paper. Below we limit ourselves to 
order-of-magnitude estimates which highlight the qualita- 
tive features of the problem. 



2.2 Radiatively inefficient accretion flow 

In this subsection we summarize the general properties of 
the radiation produced by the accretion flow. As it is men- 
tioned in the Introduction, the TeV blazars belong to the 
class of low-luminosity radio galaxies in which the accretion 
is conventi onally des c ribed in the framework of RIAF mod- 
els (see e.g. iNaravanl (|2002T ) and references therein) relevant 
for sources in which the accretion rate M is small in the 
Eddington units, 



M < Mjsdd 



0.1c 2 



: 0.2 



Mbh 



10 7 M, 



© 



M /yr 



(8) 



where MEdd is the Eddington accretion rate. 

The low radiative efficiency of the accretion flow is at- 
tributed to the low matter density which leads to inefficient 
cooling of electrons and ions in the accretion flow. As a re- 
sult, a large fraction of the released gravitational energy, in- 
stead of being dissipated in the form of radiation, is stored in 
the internal energy of matter. This energy either disappears 
under the horizon or is ejected in an outflow. 

Large internal energy of matter does not allow forma- 
tion of a geometrically thin, optically thick accretion disk. 
Instead, the accreting mat ter forms around the BH a hot op- 
tically thin gaseous torus (|Rees et al.lll982f ). In RIAF mod- 
els the accretion torus typically consists of two-temperature 
plasma with both electrons and ions being mildly relativis- 
tic. The radiation of the plasma is entirely produced by 
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electrons via three main mechanisms: synchrotron radiation, 
comptonization of the latter, and bremsstrahlung. 

The synchrotron radiation contributes to the far in- 
frared part of the spectrum. The characteristic synchrotron 



emission energy is 
heB 



^synch 



(it) 



5 x 10" 



B 



10 4 G 



1 MeV 



eV, (9) 



where B is the magnetic field in the inner part of the ac- 
cretion torus, m e is the electron mass, < 7 2 > is the mean- 
square 7-factor of electrons, T e is the electron temperature, 
and we have used that for mildly relativistic electrons 



<7e) ~ 10 



(10) 



One of the parameters of the RIAF models is the ratio j3 of 
the gas pressure p gas = nm p v 2 h /3 to the magnetic pressure 
B 2 /4n. Here n is the gas density, m p is the proton 
mass and Vth is the thermal velocity of protons which can 
be estimated as the Keplerian velocity at the corresponding 
distance, v t h(R) ~ {Rg/ R) 1 ^ 2 ■ With the use of parameter /3, 
the magnetic field is estimated as 

nl/2 



B(R) 



-^-/3 1 n(R)m p v 2 h (R) 



io 4 /r 



(11) 
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Note that we normalize the matter density in the black hole 
vicinity to the value 10 10 cm -3 , typical for the case of a 
10 7 M a black hole in the RIAF models, cf. iNaravan fc Yil 
l| 19951 ). 

If the matter density changes with the di stance as 
n(R) ~ R' 1 (in the RIAF models 1/2 < 7 sC 3/2 jLu et all 
120041 )) the magnetic field decreases as B ~ i?" (7+1)/2 , so 
that the synchrotron cooling time 
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arB 2 -y e 
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(12) 



increases with the distance as i sy nch ~ R^' +1 \ Here ot = 
0.665 x 10 -24 cm 2 is the Thomson cross-section. The syn- 
chrotron cooling is efficient only close to the BH horizon. 
Indeed, the synchrotron cooling time grows with R faster 
than the dynamical time scale given by the free-fall time, 



iff ~ \J R 3 /GMb 



(13) 



Thus, the bulk of the synchrotron emission is pro duced in 
the v icinity of the horizon. Numerical simulations (jNaravanl 
2002) confirm this qualitative result. 

The synchrotron emission in the RIAF models is 
damped due to sel f-absorption. The self-absorption coeffi- 
cient has the form l|Pacholczvkll 1970l ) 



asAU) 



where 



irhe 2 m 2 c 5 n 
2 % /3T 3 e 



2m Ve 
ZheBT} 



I(x), 



(14) 



(15) 



and the function I(x) is given explicitly in the appendix lAl 
Equating «g A to the size of the synchrotron emission region, 
which we take to be of order of the gravitational radius, one 



finds the photon energy at which the synchrotron emission 
becomes optically thin 

esA — 37e synch . (16) 

As discussed in the appendix [X] the numerical coefficient 
in this formula mildly depends on the parameters of the 
accretion flow. We omit this dependence in what follows. 

Since esa > e S ynch, the maximum of the synchrotron 
power is emitted at the energy esA- To estimate the syn- 
chrotron luminosity we approximate the spectrum by ther- 
mal radiation up to esa- This yields 



^synch 



a R 2 2T e 

^synch ^2 



(27rft) 3 



(17) 



where R 3y „ch is the size of the synchrotron emission region. 
Substituting -R syncn 
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(18) 

The part of the synchrotron spectrum above esA gives small 
contribution to the total synchrotron power because of the 
exponential cutoff in the synchrotron emission function. An 
immediate consequence of this cutoff is that the bulk of the 
synchrotron photons have energies e < esA — 0.2 eV. 

The synchrotron radiation is upscattered by inverse 
Compton (IC) process into the optical band, 



tic ~ £SA<7e) 



B 



10 4 G 



1 MeV 



eV 



(19) 



To estimate the luminosity of the accretion flow in this range 
we note that the IC volume emissivity at d istance R from 
the black hole is l|Rvbicki fc Li ghtman 2004 ) 



dLic 
dV 



-a T cn(R){~fi{R))U Bynch {R) , 



(20) 



where (7;? (7?)} is the mean-square gamma-factor of electrons 
at this distance, and J7 aync h(-R) — I/ S ynch/(47rJ? 2 c) is the en- 
ergy density of the synchrotron radiation. Integrating over 
the volume we obtain 



Lie 



— isynch O T 



n e (R){ie(R))dR . 



(21) 



Depending on the radial profiles of the electron density and 
temperature, the integral in (|21|l is saturated either close to 
the black hole horizon (for the radial profiles steeper than 
n(R){ r y 2 (R)} ~ R" 1 ) or in a region which is much larger 
than the size of the black hole. In the numerical simulations 
of RIAF models one finds the typica l size Rig of the IC 
emission region to be of order 100 R g (|Naravanll2002l ). As a 
crude estimate let us consider the case when the integral in 
(|21|) is saturated at the upper limit. Then, 



Lie ^ 0.4 x 10 3 



L. 



synch 



10 39 erg/s 
n{Ric) 



Ric 



10 8 



10 14 cm 
T e {Ri C 



1 MeV 



(22) 



erg/s 



Note that in the above estimate we took into account that 
the electron density and temperature at the distance Ric 
are lower than near the black hole. 

Finally, we consider the bremsstrahlung radiation. It 
contributes primarily into the X-ray/soft 7-ray band with 
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the total power typically comparable or lower than that 
of the synchrotron radiation. The bremsstrahlung emission 
power is proportional to n 2 , so that the bremsstrahlung lu- 
minosity from a region of the size 7? is Pbrcms ~ n(R) 2 R 3 . 
If the radial density falls down slower than n{R) ~ R~ 3 ^ 2 
(as it is the case in typical RIAF models), the bulk of the 
bremsstrahlung flux is produced at large distances. Numer- 
ical modeling of the radiatively inefficient accretion shows 
that the bremsstrahlung is produced m ainly at large dis- 
tances from the BH, Pbrcms ~ 10 5 x R g (|Naravanll2003 ). 

To summarize, the properties of RIAF important for 
the high-energy particle propagation are as follows. Matter 
distribution in a typical RIAF is characterized by a rather 
shallow radial density profile n{R) ~ P _1//2 ~R~ 3 ^ 2 , so that 
the central matter density is rather low. At the same time, 
due to the large non-dissipated mechanical energy of the ac- 
creting matter the magnetic field produced by the RIAF can 
be strong close to the central black hole {TT}. The radiation 
environment created by a RIAF has an "onion-like" struc- 
ture with the infrared synchrotron emission produced close 
to the black hole, infrared/optical IC emission produced in 
a larger region of the size ~ 100P S and X-ray /soft 7-ray 
bremsstrahlung emission produced in a region of the size 
~ 10 5 P 9 . 

To conclude this section let us make the following com- 
ment. The estimates presented above strongly depend on the 
values of the parameters of the accretion flow such as elec- 
tron temperature, density and magnetic field. These char- 
acteristics vary significantly in different RIAF models. The 
estimates presented in this section should be considered as 
indicative. A more detailed analysis of the the radiative 
background can be done using numerical simulations in each 
particular RIAF model. 

2.3 Propagation of the high-energy particles 
through the accretion flow background 

The compact particle accelerator close to the black hole pro- 
duces emission with the total power Ptot, Eq. J6}. As dis- 
cussed above, the partition of this energy between the 7-ray 
component and the accelerated matter depends on the type 
of accelerated particles, see Fig.[T] If the accelerated particles 
are electrons the power goes completely into 7-rays, while in 
the case of the proton acceleration a significant part of Ptot 
can remain in the high-energy particle beam, see Eq. ([7]). 
The primary 7-rays and protons propagate through and in- 
teract with the matter and radiation background created by 
RIAF. Interactions of high-energy particles give rise to elec- 
tromagnetic cascades that redistribute the power initially 
contained in the highest energy particles to lower energy 
bands. Since the rates of interactions of electrons and 7-rays 
are significantly different from those of protons, we consider 
these two cases separately. 

2. 3. 1 Electrons and ■y-rays 

The highest energy 7-rays produce pairs in interactions with 
the soft photon background inside the compact source. The 
cross-section of photon-photon pair-production depends on 
the center-of-mass energy of colliding photons, 

s = Ee(l - cos6)/2m 2 e , (23) 




acceleration cascade development region 

region, R ~R g R~ r 2 /^ 




10 12 10 13 10 14 10 15 10 16 10 17 D,cm 



Figure 1. Schematic representation of the two possible scenarios 
for the 7-ray emission from the compact source. Particles (elec- 
trons in the upper panel, protons in the lower panel) are initially 
accelerated in a compact region of the size of about R g , shown 
as a blue-shaded region. Particles which escape from the acceler- 
ation region initiate a cascade in the radiatively-inefficient accre- 
tion flow, at the distances R S> R g , shown as a light-blue shaded 
region. 



where E and e are the energies of the photons, and 9 is 
the collision angle. Starting from the threshold at s = 1 
the cross-section rapidly increases achieving the maximum 
(j 77 w 1.3 x 10 -25 cm 2 at s m 4, and then decreases as 
s _1 Ins. Because of relatively narrow distribution of <t 77 (s), 
7-rays interact most efficiently with the optical/infrared 
background photons of energy e ~ 1(P 7 /1 TeV) -1 eV. This 
interaction leads to absorption of the 7-rays. In order for 
the VHE emission to escape from the vicinity of the black 
hole the optical depth of this process should be less than 
one. Let us consider the contributions of the various back- 
grounds produced by RIAF to the optical depth, case by 
case. 

We start with the synchrotron background. This back- 
ground is produced in the vicinity of the black hole, P sync h ~ 
R g . Consequently, its interaction with the VHE emission 
strongly depends on the geometry of the inner part of the 
accretion disc and of the acceleration region. The situation 
is still more complicated because of the exponential depen- 
dence of th e density of the sy nchrotron background photons 
on energy |Pacholczvklll97fj| ). 

n sync h(e) = n sync h(esA)/(aO//(a:sA) > e > £ sa , (24) 

where the function I(x) is given in Eq. (|A2[) and x is defined 
in Eq. (|15p . The density n syn ch(esA) at the maximum of the 
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synchrotron spectrum can be estimated as 



2.3.2 Protons 



n-synch(eSA) = 



^synch 



^synch^AC 



(25) 



10 1 



synch 



10 40 erg/s 



ynch 



10 12 cm 



T esA 1 - 1 
Lo.2 cvJ C 



Using Eq. (|24|l one finds that because of the sharp cut- 
off in the synchrotron spectrum above the energy esA ~ 
0.2 eV, 7-rays with energy E 1 ~ f TeV escape through 
the synchrotron background, r 77 (i5 7 = ITeV) ~ f, if 
J^synch ~ f0 12 cm and L aync ^ < f0 40 erg/s. On the other 
hand, our crude estimates indicate that the spectrum of 
the escaping 7-rays should be sharply cut-off at the en- 
ergy E^ — 5[e S A/0.2 eV]" 1 TeV, since t 77 rapidly grows 
to r 77 2> 1 at this energy. One should note, however, that 
the details of the behaviour of the spectrum close to the cut- 
off strongly depend on the details of the spatial distribution 
of the synchrotron emission. Since the threshold of the pair 
production, s — 1 (see Eq. (|23p 1. depends on the angle 6 
between the 7-ray and synchrotron photon velocities, the 
anisotropy of the synchrotron emission can result in a shift 
of the cut-off in the 7-ray spectrum toward higher energies. 

VHE 7-rays produced outside of the synchrotron emis- 
sion region, or those which escape from it, pass through the 
IC background with the photon density: 



ftic 



Lie 



47ri? 2 c e 



(26) 



1.7 x 10 1 



10 39 erg/s 



Ric 



10 1 



It 



cVJ 



The optical depth of this background for the TeV 7-rays can 
be estimated as 



T 77 (Ery) 



■■ <7 77 nic-Ric 



2.2 



10 39 erg/s 



Rk 



101 



En, 



1 TeV 



(27) 



Thus, the inner part of the accretion flow is transparent for 
the TeV 7-rays if 



Lie < 0.5 x 10 3 



Ric 



10 14 cm 



1 TeV 



erg/s . (28) 



Finally, the background of the bremsstrahlung radiation 
does not affect propagation of the high-energy 7-rays. This 
is due to two reasons: the low density of the corresponding 
photons and their relatively weak interaction with the VHE 
7-rays because of the l/s suppression of the photon-photon 
cross-section. 

Thus, in the case when the luminosity of the accretion 
flow is as low as L acc < 10 40 erg/s, the VHE 7-ray emission 
can originate directly from the vicinity of the black hole, ft 
can be e.g. the direct synchrotron/curvature emission which 
accompanies the acceleration process. This possibility is il- 
lustrated schematically in Fig. 1 (upper panel). Otherwise, 
if the synchrotron luminosity of RIAF is significantly above 
10 40 erg/s and/or the IC luminosity exceeds the limit (I28|) . 
the VHE 7-ray emission must be produced outside the cor- 
responding radiation regions. This is possible if the primary 
accelerated particles are protons. 



Protons with energies above the threshold 



E t h = 



m 7r (m 7r + 2mp) 



4c 



8 x 10 1 



[l cvi 



eV (29) 



loose energy in the interactions with the infrared photons 
of energy e via pion production. Near the threshold, the 
cross-section is dominated by the single pion production res- 
onance <t P7 ~ 6 x 10 -28 cm 2 ; in this regime proton gives 
only 20% of its energy to the pion. At higher proton ener- 
gies the photo-pion production cross-section decreases down 



to o r . 



10" 



cm , while the proton energy loss in every 



interaction increases up to 50%. Thus in both cases protons 
lose most of their energy at similar distance. 

For high synchrotron background, L sync h ~ 10 42 erg/s, 
protons with energies E p > 10 18 eV cannot escape from the 
acceleration region, 



synch 



PI 



(E p > 10 18 eV) = (Jp-y fisynch R; 



synch 



L 



synch 



10 42 erg/s 



R 



synch 



10 



12 



cm 



[o.2 cvi ' 



(30) 



On the other hand, protons with energy E p < 10 eV inter- 
act only with exponential tail of the synchrotron emission 
and escape from this region. 

During the propagation through the inverse Compton 
emission region at the distance scales R ~ 100i? 9 , protons 
interact with the inverse Compton photon background, so 
that the optical depth is 



ic _ 



0.2 



P7 nic-Ric 
Lie 



10 41 erg/s 



Ric 



10 14 cm 



[l ev] 



(31) 



The power of the absorbed proton flux is converted into the 
products of pion decays: neutrinos, 7-rays and electrons of 
the energies ~ 0.1-Ep. Obviously, the neutrinos freely escape 
from the production region. 

Naively, from the above discussion of propagation of 7- 
rays, one would expect that 7-rays are not able to escape 
from such a dense IC background. However, this is incor- 
rect. The point is that the energies of the 7-rays produced 
in the neutral pion decays are much above the energy cor- 
responding to the maximum of the pair production rate on 
the IC background. For example, for the 10 17 eV 7-rays the 
cross-section of the pair production on the IC background 
is s/lns ~ 10 4 times smaller than the peak value, so that 
it is of the order of p-y interaction cross-section. Thus, the 
mean free path of the secondary 7-rays is comparable to the 
mean free path of the primary protons and is of order of the 
size of the IC radiation region. These 7-rays can give rise 
to electromagnetic cascades redistributing energy into TeV 
7-rays. The draw back of this mechanism of TeV emission is 
its low efficiency (jNeronov fc Scmikoz 2003) due to the fact 
that the TeV 7-rays which are observed at infinity can be 
produced only in the surface layer of the IC radiation re- 
gion, where the optical depth for the TeV radiation drops to 
the values of order one. On the other hand, this mechanism 
should result in strong multi-GeV emission which can escape 
through the whole inverse Compton radiation region without 
significant absorption. This implies that if this mechanism 
is indeed responsible for the TeV emission in (some of) the 
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TeV blazars, it should lead to a strong signal in the multi- 
GeV band accessible for the Fermi (GLAST) satellitsfl. 

Apart from the interactions with soft radiation back- 
ground, protons can also interact with the matter in the jet. 
The optical depth of protons with respect to this process is 
estimated as 



; a pp n p (R)R . 



(32) 



where n p (R) is the matter density in the jet at the distance 
R. Taking the pp interaction cross section at energies E p ~ 



10 eV to be equal to a p 



10" 



0.1 



n p (R) 



10 s cm 



cm we obtain 



R 



10 16 cm 



(33) 



Particle multiplicity in pp collisions at E v ~ 10 18 eV is about 
N ~ 100 — 200 (see e.g. (|Heiselberdl200ll ) for a review) . Thus 
every collision of a 10 18 eV proton with a background proton 
produces ~ 100 photons with E 1 = 5 x 10 15 eV from 7r°, 
~ 100 neutrinos with the same energy from 7T decays and 
~ 100 electrons and positrons. All neutrinos escape from the 
interaction region, while electrons, positrons and photons 
give rise to an electromagnetic cascade. As we discuss in the 
next section, this provides a link with the standard picture of 
VHE emission by relativistic blobs in the jet. The limitation 
of this mechanism is that it requires rather large densities 
of matter in the jet, see Eq. (|33|l . 

To summarize, the qualitative analysis of this subsec- 
tion shows that production of VHE 7-rays can, in principle, 
be possible within the compact source model with proton 
acceleration, see Fig. 1, lower panel. Still, the details of the 
mechanism of conversion of the proton flux power into the 
power of TeV 7-ray emission are yet to be worked out. 



2.4 Particle cascade in the accretion flow: a link 
to the "relativistic blob in the jet" picture 

The 77, P7 and pp interactions taking place during the prop- 
agation of the high-energy particles through the RIAF en- 
vironment give rise to electromagnetic cascades, which re- 
distribute the power initially contained in the highest en- 
ergy particles to the lower energy bands. This is expected 
to result in the broad-band (radio-to-7-ray) electromagnetic 
emission from the cascade. 

To large extent, this emission can be described by 
the conventional synchrotron - self Compton and/or syn- 
chrotron - external Compton models based on the picture 
of relativistically moving blobs of plasma in the jet. Indeed, 
as already mentioned, the primary flow of high-energy par- 
ticles, accelerated by the compact source, is expected to 
be highly anisotropic, with the direction set by the mag- 
netic field in the acceleration region. The cascades produce 
a stream of relativistic particles with velocities scattered 
within some angle £ around the direction of the primary 
flow. From the kinematical point of view such a stream is 
nothing else than a blob of plasma moving with the bulk 
Lorentz factor 



r ~ - 



(34) 



2 The official website of the Fermi collaboration is 
http: / / fermi.gsfc.nasa.gov/ 



The angle Q is determined by the dynamics of the cascade. 
An important role in this dynamics is played by the value 
and configuration of the magnetic field. At the late stages 
of the cascade development, when the density of the par- 
ticles in the cascade is large, the problem should be solved 
self-consistently taking into account the back-reaction of the 
plasma in the cascade on the magnetic field. This observa- 
tion makes the link between the cascade and the blob pic- 
tures not only kinematical but also dynamical: in the stan- 
dard approach the magnetic field is also determined self- 
consistently by the dynamics of the blob itself. 

Existence of this link enables to use in the analysis of 
TeV 7-ray flares within the framework of the compact source 
models many results from the standard approach. In partic- 
ular, the transparency of the blob for the TeV radiation in 
the Jul y 2006 large flare of P KS 2155-304 implies the con- 
straint jBegelman et al.ll2008h T > 50 and hence, £ < 1.1°. 

On the other hand, there are several important differ- 
ences between the cascade and a generic relativistic blob 
model which potentially enable to distinguish the two mod- 
els observationally. First, in the cascade setup, the relativis- 
tic blobs are formed in a compact region close to the black 
hole. In this way this setup naturally incorporates the ob- 
served fast variability and is not in a direct conflict with 
the low bulk Lorenz factors observed at the parsec-scale dis- 
tances. Indeed, the bulk motion of the secondary particles 
in the cascade can decelerate at parsec-scale distances either 
because of development of intrinsic instabilities or because 
of interaction with the interstellar medium. To study this 
possibility, a dynamical numerical calculation of the prop- 
agation of the high-energy particle cascade through the ac- 
cretion flow and through the interstellar medium is needed. 
Second, the cascade transfers power from the higher-energy 
particles to the lower-energy ones: in this sense it is a "top- 
down" scenario of formation of the spectrum of emitting 
electrons. It would be interesting to understand if the latter 
property can explain the existence of the low-energy cut- 
offs in the spectra o f electrons in the s ynchrotron - inverse 
Compton scenarios (Krawczvnski 2007). Finally, as already 
mentioned, in the compact source model, the emission from 
the blob may be superimposed onto the 7-ray component 
originating directly in the acceleration region near the black 
hole. 



3 TIMING OF THE SUPERMASSIVE BLACK 
HOLE IN PKS 2155-304 WITH TEV 7-RAYS 

Within the compact source model the spectral and variabil- 
ity properties of the VHE signal can be interpreted in a 
qualitatively different way (as compared to the model of 
VHE emission from the parsec-scale distances). They can 
be used to "probe" the physical environment in the central 
engine: the matter density, the radia tive efficiency of the 
accretio n flow, the magnetic field (cf. iRieger fc Aharonianl 
((2008bJ)). Moreover, in this case the VHE observations pro- 
vide a completely new tool for the study of the black hole 
physics. Namely, the 7-ray timing data can be used to derive 
constraints on the parameters of the black hole itself. Below 
we explore such a possibility on a particular example of the 
July 2006 la rge flare of PKS 2155- 304 (redshift z = 0.116) 
reported by l|Aharonian et al1l2007l ). 
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3.1 The relevant time scales 

Timing properties of emission produced in the vicinity of 
black hole horizon can be characterized by several funda- 
mental time scales. First, the minimal possible variability 
time scale is determined by the requirement of causal con- 
nection of the emission region and is given by the light cross- 
ing time of the black hole horizon, 



(35) 



Rg 




tic = 2 (ll g + ytRI-a 2 ^ /c 

10 2 [Af BH /lO 7 A/ ] s, 
2 x 10 2 [A/ BH /1O 7 M ] s, 

where R g is defined in Eq. ([2|; the parameter a is related 
to the angular momentum Jbh of the black hole as a = 
Jbh/Mbhc 2 and lies in the range ^ a ^ R g . 

If the 7-ray emission is produced close to the black hole, 
variability at the characteristic time scale of rotation around 
the black hole is expected on general grounds, unless the 
entire system (the accretion flow with an embedded parti- 
cle acceleration region) is perfectly axially symmetric. The 
intensity of the modulation of the signal can depend on vari- 
ous parameters, such as the inclination angle of the observer 
with respect to the BH rotation axis, distance of the emis- 
sion region from the horizon, etc. A perfect axial symmetry, 
which would wash out the modulation of the signal with the 
period of rotation around the black hole, can be expected 
in a stationary, quiet state of the source. On the contrary, 
a bright flare is, most probably, related to rapid change of 
the system parameters (for example, inspiralling of a denser 
clump of matter into the black hole) which should lead to a 
significant disturbance of the axially symmetry. 

Close to the horizon of a spinning black hole, the ac- 
creting matter rotate s around the black hole with the period 
ijBardeen et alJll972h 

P(r)=2« r tj° , (36) 

CKg 

where r is the radius of the orbit. The + (— ) sign corre- 
sponds to the prograde (retrograde) orbit. It is known that 
stable circular orbits exist only down to certain distance r ms 
from the BH. The period of rotation along the last prograde 
stable orbit at the distance r ms is 



P(r ms ) 



inRg/c- 630 yffife] s, 
!2V6nR g /c ~ 4600 



a — R g 
a = 



(37) 

The disturbance of the axial symmetry of the accretion flow 
is expected to result in the modulation of physical conditions 
in the AGN central engine with a period given by Eq. (I37p . 
Since the properties of the 7-ray emission from the central 
engine (from the base of the jet) are determined by the phys- 
ical conditions in the central engine, the modulation of these 
conditions should result in the modulation of the 7-ray sig- 
nal with the same period. 



3.2 Fitting the lightcurve of PKS 2155-304 

The VHE 7 -ray lightcurve of the flare of PKS 2155-304, 
reported in ()Aharonian et all 12007m . consists of several pro- 
nounced sub- flares (see Fig. [2} . At least three characteristic 
time scales can be found in a straightforward way from the 




30 40 50 60 70 80 90 100 110 120 130 
t [min] 

Figure 2. The lightcurve of PKS 2155-304 fitted with a sequence 
of sub-flares with identical profiles, but with different normaliza- 
tions (see Table [TJ. 



analysis of the lightcurve: the rise and decay times of the in- 
dividual sub-flares, irisc , tdecay , and the period of recurrence 
of the sub-flares, T. In order to estimate the average val- 
ues of these parameters we have fitted the VHE lightcurve 
assuming that the time profiles of the individual sub-flares 
are characterized by the same rise and decay times, and 
differ only in the arrival times and the amplitude^- The ap- 
proximation of identical i r ise> tdecay appears naturally in the 
context of the compact source models where these parame- 
ters are determined by the intrinsic time scales of the central 
engine, such as e.g. the light crossing time of the black hole 
and the cascade development time. On the other hand, it 
is clearly a strong idealization: the complicated dynamics 
of the particle acceleration and the cascade development is 
expected to introduce a scatter in the characteristics of the 
individual sub-flares. Still, we stick to this approximation as 
a natural first step. 

The individual sub-flares are modeled with the profile 

T U\ J Nk ex P[(* ~ *max,fc) Arise] , t < t 

max, re /oo\ 
[Nk exp[— (t — £max,fc) Adccay] , £ > tmax.fc 

where the time i ma x,fc corresponds to the maximum intensity 
and Nk is the amplitude of the fc-th sub-flare. Note that we 
define i r isc (idecay) as the time in which the signal increases 
(decreases) by a factor e. We fit the lightcurve with the sum 
of several sub-flares (|38[) and a constant signal. The result of 
the fit of the overall lightcurve with such a model is shown 
in Fig. [5] The rise and decay times inferred from the fit are 



t lisc = (2.5 ± 0.2) x 10 2 s, 
tdocay = (4.9 ± 0.5) x 10 2 



(39) 
(40) 



The arrival times t mlXj t and normalizations Nk of the sub- 
flares are summarized in the Table [T] 

Initially, we fitted the lightcurves with a set of five sub- 
flares which correspond to the five pronounced peaks clearly 
visible in the data. The quality of this fit was rather low 

3 This model is diff erent from the model considered by 
jAharonian et al]l200^ where the rise and decay times were al- 
lowed to vary among individual sub-flares. 



© 0000 RAS, MNRAS 000. IM 



Measuring parameters of AGN central engines with very high energy ^-ray flares 9 



k 


£max,fc[min] 


A r fc [10 _9 cm- 2 s- 1 ] 


I 




2.8 ± 0.3 


2 


46.2 ±0.8 


0.78 ± 0.24 


3 


59.1 ±0.2 


3.0 ±0.17 


4 


69.8 ±0.3 


2.9 ±0.2 


5 


79.0 ±0.4 


1.8 ±0.2 


6 


89.9 ±0.3 


2.5 ±0.2 


9 


120.7 ± 1.0 


0.30 ±0.11 



Table 1. List of parameters of the sub-flares. 

(X 2 = 106 for 70 de grees of freedom). We have found that 
the addition of one more sub-flare (sub-flare number 2 in 
Table [TJ significantly improves the quality of the fit (\ 2 = 
86.3/68 d.o.f.; the F-test gives a chance probability of the fit 
improvement at the level of 0.1%). Finally, we have found 
that the fit can be further improved by addition of one more 
sub-flare near the end of the overall lightcurve, the last sub- 
flare in the Table [T] The quality of the fit with seven sub- 
flares is x 2 = 75.8/66 d.o.f. According to the F-test the 
probability that the latter fit improvement is achieved by 
chance is 1.4 %. 

The arrival times of the six bright sub-flares follow an 
approximate linear law 



imax,fc — to ± k • T 



(41) 



see Fig. [3] This is reminiscent of quasi-periodic oscilla- 
tions observed in the X-ray band in the X-ray binaries 
( van der Klislk oOO) and in the infrared band in the Galac- 
tic Center (jGenzel et al]|2003l ). The last weak sub-flare also 
falls on the linear dependence (|41|l if assigned the number 
k — 9. However, we do not use this sub-flare in the following 
analysis because its presence is not strongly required by the 
lightcurve fit. 

Fitting the set of arrival times t max .fc of the six bright 
sub-flares (numbered by 1-6 in the Table [1} with the linear 
law (|41[) one finds the recurrence period of the sub-flares, 



T= (5.93 ±0.15) x 10 2 s. 



(42) 



This recurrence period coincides wi th the characteris t ic var i- 
ability scale ~ 600 s mentioned bv lAharonian et alJ (|2007T I. 

We have tested if the apparent quasi-periodicity of the 
signal can appear in the data by chance. To do this we have 
fixed the arrival times of the first and last of the bright 
sub-flares (sub-flares 1 and 6 in the Table [TJ. Then we have 
allowed the arrival times of the 4 intermediate sub-flares 
to be distributed randomly between t max ,i and t max ,6- The 
arrival times i maXj fc, k = 1, . . . , 6 are assigned the errors from 
the Table [T] We have calculated the probability that fitting 
the arrival times of the six first sub-flares by the function 
(|4ip would result in a y 2 ^ Xo> where Xo = 41 is the x 2 



of the fit of the real datqj. This probability turns out to be 
1.3 • 10 -3 . We also tried another error assignment. Namely, 
all the arrival times were assigned the error St — 1.5 min, 
which corresponds to the intrinsic scatter of the arrival times 

4 Note that Xo is large. This means that the oversimplified model 
of strictly identical sub-flares with strictly periodic arrival times 
is actually excluded by the data. This is not surprising: as already 
mentioned, on physical grounds one expects deviations from this 
idealization. 




Figure 3. The arrival times of the sub-flares of the PKS2155-304 
lightcurve as a function of the sub-flare number. The data are 
fitted by the straight line t max k = *0 + & ■ T. 



of the real sub-flares around the linear law (|4ip . In this case 
the best linear fit corresponds to xi = 5.43 for 4 degrees 
of freedom. The chance probability to obtain x 2 ^ 5.43 in 
the simulated data sets with random arrival times of the 
intermediate sub-flares turns out to be 6 ■ 10~ 3 in this case. 
We have checked that our conclusions do not depend on the 
particular choice (|38[) of the sub-flare profile by considering 
other possible choices of Ik(t)- 



3.3 Constraints on the black hole parameters 



The rise and decay times (|39[) . (|40p as well as the recur- 
rence time (|42p are determined by the physics of the 7-ray 
emission. Within the compact source scenario described in 
Sec. [2] the minimal possible rise time is given by the light 
crossing time of the black hole, Eq. (|35[1 . Requiring 



trisc ^ (1 + Z)t\ c 



(43) 



one finds a range of parameters Mbh , a, excluded by the ob- 
servations. This range corresponds to the dark shaded region 
in Fig. H 

This should be compared with the constraints obtained 
by other methods. In the literature there are two estimates 
of the mass of the central black hole in P KS 2155-304. 

The first one (jAharonian et al.ll2007h is obtained from 
the relation between the masses of ce ntral black ho l es an d 
luminosities of host galaxy bulges ()Bettoni et al . 2003); 
it gives Mbh ~ 10 9 Mp. As it was already noted in 
|Aharonian et al]|2007t ). this estimate is in conflict with the 
constraint Mbh ^ 2 x 1O 7 M0 derived from the 7-ray vari- 
ability properties. We point out, however, that there are sev- 
eral uncertainties in the estimate based on the Mbh — ibuigo 
relation. It is obtained by extrapolating the empirical re- 
lation observed in a local sample of normal galaxies to 
the case of TeV blazars. The reliability of this extrapola- 
tion was never investigated. Besides, the relation itself has 
a large intrinsic scatter (order-of-magnitude deviations are 
present). Finally, there is no consensus in the literature 
abou t the luminosity of th e host galaxy in PKS 2155-304 
(see l)Aharonian et al.ll2007h and references th erein). 



The second estimate jzhang et al.l 120051 ) is based on 



© 0000 RAS, MNRAS 000, [THE] 



10 Neronov et al. 




10° - 



1 

ac 3 /GM 

Figure 4. The parameters of the supermassive black hole in PKS 
2155-304 inferred from the VHE 7-ray lightcurve. The shaded 
region shows the range of JIbh, o, excluded by the requirement 
1143 1 . Thick blue curve shows the relation between Mbh and a 
obtained by identifying the period of the sub-flare recurrence with 
the rotation period over the last stable orbit. The horizontal line 
shows the lower bound on t he black hole mas s implied by the 
X-ray variability analysis of llZhang et 3,1.11200511 . For a > this 
bound should be taken with caution as its dependence on a has 
not been explored. 



the study of X-ray var iability properties of P KS 2155-304: 
it uses the method of ijNikolaiuk et a.1.1 l2004h to relate the 
mass of black hole to the excess variability a„ xs at a cer- 
tain frequency scale. This method yields the bound Mbh 1? 
8.1 x 1O 6 M0, shown by a dashed horizontal line in Fig. [4] 
The latter estimate is compatible with the constraint derived 
from the TeV variability of the source. One should note, 
however, that the estimate of the black hole mass based on 
the X-ray variability properties suffers from the same uncer- 
tainty as the estimate based on the Mbh — ibuigc relation: 
it is originally derived for a sample of nearby non-blazar 
AGNs. Its applicability to the sample of the TeV blazars 
was never tested. 

The tight constraint on the black hole mass, derived 
from the TeV variability time scale, provides a possibility 
of strong observational test of the compact source model. 
Indeed, a precise determination of the black hole mass by 
another method would be able to falsify or confirm the com- 
pact source model. 

The indication of the quasi-periodicity of the sub-flare 
arrival times suggests to associate the period T of the sub- 
flare recurrence with the (minimal possible) period of rota- 
tion around the black hole, 

T = (l + *)P(r mB ). (44) 

This gives a relation between a and Mbh shown by the thick 
blue curve in Fig. [4] If combined with the constraint on the 
bla ck hole mass deriv ed from the X-ray variability analy- 
sis (|Zhang et al.ll2005T l. this relation implies that the black 



hole is rotating almost at the maximal rate, a ~ R g . How- 
ever, we remind that the X-ray constraint on the black hole 
mass should be taken with caution. In particular, to the best 
of our knowledge, the dependence of this phenomenological 
constraint on the black hole spin a has not been explored. 



4 CONCLUSIONS 

In this paper we have proposed that the recently observed 
fast variability of the VHE emission from blazars can be 
naturally accommodated within the framework of "compact 
source" model. In this model particles responsible for the 
observed VHE emission are accelerated close to the cen- 
tral supermassive black hole, rather than at large distances 
downstream the AGN jet. We have analyzed the problem of 
escape of the VHE 7-rays from the vicinity of the central 
engine and demonstrated that the region around the central 
engine is transparent for TeV 7-rays if the accretion flow 
in the TeV blazars is radiatively inefficient. If the luminos- 
ity of the accretion flow is as low as L acc < 10 40 erg/s, the 
TeV 7-ray emission can come directly from the immediate 
neighborhood of the central black hole. Alternatively, for 
brighter accretion flows, the TeV 7-rays may be produced 
at some distance from the black hole in a proton-initiated 
electromagnetic cascade developing in the accretion flow. 

The possibility that the properties of the VHE 7-ray 
emission are directly linked to the properties of the cen- 
tral engine of the AGN, if confirmed by future observations, 
provides a new tool to study the physical conditions in the 
direct vicinity of the supermassive black hole. In particular, 
the VHE signal can be used to constrain the parameters of 
the accretion flow and of the black hole itself, such as its 
mass and spin. 

We demonstrated this possibility on the example of the 
bright TeV flare of the blazar PKS 2155-304. Within the pro- 
posed scenario, the characteristic time scales, found in the 
timing analysis of this flare, are directly related to the pa- 
rameters of the supermassive black hole. The minimal vari- 
ability time scale of the signal is identified with the black 
hole light-crossing time. This sets the bound on the black 
hole mass and its rotation moment shown in Fig. [3] We also 
observed that the signal exhibits quasi-periodic oscillations. 
Identifying the recurrence time of these oscillations with the 
period of rotation around the black hole we obtained a rela- 
tion between the black hole mass and its rotation moment. 

A detailed modeling based of the framework proposed in 
this paper should involve calculation of particle acceleration 
and propagation in the vicinity of the black hole through 
the environment created by the accretion flow. This can be 
done assuming particular (numerical) models of RIAF and 
particle acceleration. We leave this for future work. 
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APPENDIX A: SYNCHROTRON 
SELF-ABSORPTION IN THE ACCRETION 
FLOW 

The function I(x) entering into the expression (|14[l for 
the synchrotron self-absorption coefficient has the form 
(|Pacholczvklll970D . 



1 f c 



z F{x/z 2 )dz , 



(Al) 



where F(x) = x J°° K 5 / 3 (z)dz. It is straightforward to ob- 
tain the asymptotics of I(x), 



/2/37rexp (- 3(x/4) 



l/3\ 



X > 1 



(A2) 



The synchrotron photons are self-absorbed below a cer- 
tain energy esa- The latter is estimated from the condition 
that the optical depth for the synchrotron emission with 
energy esA is equal to unity, 



«SA(eSA)-Rsynch = 1 



(A3) 



where -R syn ch is the size of the synchrotron emission region. 
Taking i? sync h fa R g one obtains the following equation for 
the variabl^EI ^sa corresponding to the self-absorption en- 
ergy e S A, 

1.894a +lna;sA = 17.37 



+ ln 



n 



10 10 cm- 
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_WM q _ 




_10 4 G_ 




1 MeV_ 



(A4) 



When the logarithm in the second line vanishes, the solution 
to this equation is 



x S a = 247 



(A5) 



This corresponds to the value (|16jl of the self-absorption 
energy. Note that the r.h.s. of the equation (|A4|) logarith- 
mically depends on the parameters of the accretion flow im- 
plying that the variable :tsa is only mildly sensitive to these 
parameters. Namely, the value (1A5[) is multiplied by a factor 
ranging from 0.6 to 1.5 when the combination of the param- 
eters entering the logarithm in the second line of Eq. ()A4|) 
varies from 0.1 to 10. 
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